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We study the QED bound-state problem in a light-front hamiltonian approach. Starting with a bare 
cutoff QED Hamiltonian, H B , with matrix elements between free states of drastically different energies 
removed, we perform a similarity transformation that removes the matrix elements between free states 

. 

with energy differences between the bare cutoff, A, and effective cutoff, A (A < A). This generates effective 

<N . 

interactions in the renormalized Hamiltonian, H„. These effective interactions are derived to order a in 

CO 

this work, with a -C 1. H R is renormalized by requiring it to satisfy coupling coherence. A nonrelativistic 

m 

limit of the theory is taken, and the resulting Hamiltonian is studied using bound-state perturbation 

in 
o 

\o 

0\ 



theory (BSPT). The effective cutoff, A 2 , is fixed, and the limit, < — m 2 a 2 <C A 2 <C m 2 a — ► oo, is 
taken. This upper bound on A 2 places the effects of low-energy (energy transfer below A) emission in the 
effective interactions in the |ee) sector. This lower bound on A 2 insures that the nonperturbative scale of 
interest is not removed by the similarity transformation. As an explicit example of the general formalism 



introduced, we show that the Hamiltonian renormalized to 0(a) reproduces the exact spectrum of spin 



, — i 

x. 

splittings, with degeneracies dictated by rotational symmetry, for the ground state through 0(a 4 ). The 
entire calculation is performed analytically, and gives the well known singlet-triplet ground state spin 
splitting of positronium, -^a 2 Ryd. We discuss remaining corrections other than the spin splittings and 
how they can be treated in calculating the spectrum with higher precision. 
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I. INTRODUCTION 



There is much effort being put into solving for the hadronic spectrum from first principles of QCD in 
(3 + l)-dimensions using a light-front similarity hamiltonian approach [1-5]. However, low-energy QCD is 
challenging, and a realistic analytical calculation may be impossible. There is a need for exact analytical 
calculations that test and illustrate the approach. This paper is one such calculation. The calculation 
in this paper was discussed in Ref. ||, where the leading order calculation was completed. We start 
from the canonical QED 3+1 Hamiltonian, and set up the general formalism for deriving the renormalized 
Hamiltonian, H R , to some prescribed order in a (with q«1), and then using BSPT to compute the low 
energy observables of interest to some prescribed order in a. 

The general formulae are then applied explicitly to the spin splittings in positronium. Renormalization 
is carried out to order a, and a nonrelativistic limit of the theory is taken. The limit , 

< m 2 a 2 < A 2 < m 2 a ► oo , (1) 

is taken, and with BSPT, a calculation of the singlet-triplet ground state spin splitting of positronium 
to order a is carried out. The entire calculation is performed analytically, and the well known result, 
^a 2 Ryd, is obtained. Our results apply to all physical positronium states with arbitrary momentum V . 

Kaluza and Pirner have calculated the singlet-triplet ground state spin splitting of positronium (ne- 
glecting annihilation channel contributions) to order a 4 in a light-front hamiltonian approach || , and they 
obtained correct results numerically, but were forced to make ad hoc assumptions because their Hamilto- 
nian depended on the full eigenvalue of the problem. We avoid these assumptions in our approach, and 
perform the calculation analytically. 

II. THE SIMILARITY HAMILTONIAN APPROACH 

The starting point in the similarity hamiltonian approach is a bare cutoff continuum Hamiltonian 
of physical interest, H B , with energy widths^ restricted to be below the bare cutoff, A. A similarity 

2 The "energy difference" between the free states in a matrix element of a Hamiltonian is defined to be its "energy width." 
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transformation (unitary here) is defined that acts on H B and restricts the energy widths in the final 
Hamiltonian, H R , to be below the effective cutoff, All the energy changes between A and A are 
"integrated out" and replaced by effective interactions in H R . The initial Hamiltonian must then be 
adjusted so that H R satisfies coupling coherence Q, which produces a renormalized Hamiltonian order by 
order in the running couplings, e\ and m\. Coupling coherence can be realized by requiring a Hamiltonian 
at one scale to equal a Hamiltonian at a new scale after changing the explicit scale dependence in the 
Hamiltonian and the implicit scale dependence in a finite number of independent running couplings. In 
addition, all the dependent couplings are required to vanish when the independent marginal couplings 
are taken to zero. 

The second step of this similarity hamiltonian approach is the diagonalization of H R . First, for this 
QED calculation, a nonrclativistic limit of H R is taken. This is reasonable because the bound-state 
electron momenta (equivalent formulae hold for the positron) satisfy: 

Pel^™ = I + 0(a), (2) 

Pe ' ectro " = Oia) , (3) 

TO 

a < 1 , (4) 

where to is the renormalized electron mass and V + is the total longitudinal momentum of positronium.^] 
Next, the Hamiltonian is divided into a nonperturbative and perturbative part, 

H R = H + (H R -Ho)=H +V . (5) 

This is a standard trick for hamiltonian problems, and will work best if the lowest order spectrum of H R 
is well approximated by the spectrum of Ti . Phcnomcnological input can be used to determine Ti. . The 
main point is that if the spectrum of Tl differs too much from the lowest order spectrum of H R , the 
subsequent BSPT will not converge rapidly (if at all). 

3 In the initial setup of the similarity transformation, A and A will be used as a shorthand for ^jqr and respectively, 

where P+ is the total longitudinal momentum of the physical state of interest, and A 2 and A 2 have dimension (mass) 2 . 
4 The reader unaccustomed to light-front coordinates should consult Appendix A. 
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We proceed to estimate a lower bound on A where the nonperturbative bound state effects of H\ 
enter. This lower bound naturally appears at the physical binding energy of interest: 

A 2 |M^-(2m) 2 | 

r+ » J r+ L ' 

where Mn is the mass of the physical state of interest. 5 Given this lower bound we next proceed with 
an estimate of an upper bound on A that will allow us to obtain our approximate spectrum entirely 
within the |ee) sector. The fact that this upper bound arises at all is one of the utilities of the similarity 
hamiltonian approach. This upper bound is 

A 2 

p+ < Won > ( ? ) 

where Q~ hoton is the dominant energy of emitted or absorbed photons. 

Given this range of A that allows the mass eigenstates of the theory to be dominated by the few-body 
sector of interest, and at the same time, does not remove the nonperturbative bound state physics of 
interest, we proceed with the BSPT and calculate all relativistic corrections to the zeroth order spectrum 
to some specified order in a. The utility of this paper is that the relativistic corrections to order a 4 in 
the ground state spin splittings can be calculated analytically. In the following subsections we describe 
the similarity hamiltonian approach more explicitly. 

A. Step one: derivation of H R 

A self-contained discussion of the derivation of the renormalized Hamiltonian, H R , will now be given. 
The similarity transformation we use to derive our effective interactions was developed by Glazek and 
Wilson Q , and separately by Wegner || . An early application of this approach was a weak-coupling 
treatment of QCD. 

The approach starts with a definition of the bare Hamiltonian: 

H B = h + v A , (8) 

5 M^ = (2m + Bjv) 2 , where N labels all the quantum numbers and — -Bjv is the binding energy of the physical state of 
interest . 



V A = A^A . (9) 

v A = v can + Sv A , (10) 

H can = h + V can , (11) 

where h is the free Hamiltonian, H can is the canonical Hamiltonian, f A is a regulating function, and 
Sv A are counterterms defined through the process of renormalization. The canonical Hamiltonian, H can , 
is written in terms of rcnormalized parameters and will be specified at the end of this section. The 
counterterms, Sv A , are fixed by coupling coherence. Coupling coherence will be explained further below. 
The free Hamiltonian, h, is given by: 

h = I E { ( b t(P) b *(P) + dt(p)d s (p)) p "t + a l(P) a °(P) (j^j } ■ (") 

h\i)=s i \i) , ^2\i)(i\ = l, (13) 

i 

where the sum over i implies a sum over all Fock sectors and spins, and integrations over all momenta 
in the respective free states. We use the shorthand J p = J d - ' i^3p+ P ' ■ 171 1S the renormalized fermion 
mass. 

The regulating function, J A , is defined to act in the following way: 

(*I/a^aIJ> = /Ay(^ A |i) =f Aij v Aij , (14) 
f Aij = tf(A-|Ay|) , Ay = £i - ej . (15) 

Note that this choice of a step function is not necessary and can lead to pathologies, however it is useful 
for doing analytical calculations. 

Next, a similarity transformation is defined that acts on H B and restricts the energy widths in the 
renormalized Hamiltonian, H R1 to be below the effective cutoff, A. This transformation allows recursion 
relationships to be set up for H R , which can be written in the following general form: 

H R = h + v Xl (16) 
v x = f\v\ , (17) 

_ _(1) _(2) 

vx = v x +v x +■■■ , (18) 



where the superscripts imply the respective order in v can . 

Now, starting with the above bare Hamiltonian, we will describe this procedure more explicitly. The 
similarity transformation is defined to act on a bare cutoff continuum Hamiltonian, H B , in the following 
way: 

H R = S(X,A)H B S^(X,A) , (19) 
S(X,A)S^(X,A) ee 5 t (A,A)5(A,A) ee 1 . (20) 

This transformation is unitary so H B and H R have the same spectrum: 

H B \* B ) = E\9 B ), (21) 
S(X,A)H S\X,A) S(X,A)\* B ) = ES(\ t A)\Hf B ) . (22) 

Therefore, E is independent of the effective cutoff, A, if an exact transformation is made. E is also 
independent of the bare cutoff, A, after the Hamiltonian is renormalizcd. 

To put the equations in a differential framework, note that Eq. ( [l9| ) is equivalent to the following 
equation 

dH R 
dX 

with 



[H R ,T X ] , (23) 



S(X,A) ee TexpH Tx<dX^J 



(24) 



where 'T' orders operators from left to right in order of increasing energy scale, A'. Eq. (E3) is a first 
order differential equation, thus one boundary condition must be specified to obtain its solution. This 
boundary condition is the bare Hamiltonian: H R \ X ^ A ee H b . H b is determined by coupling coherence. 
Now we must specify T\, the anti-hermitian (rj = — T\) generator of energy width transformations. To 
define T\ note that it is enough to specify how v\ and h change with the energy scale A. This is seen by 
writing out Eq. ( |23| ) more explicitly using Eq. jl^): 

d^ + lX ihlJx) = {h,T x } + [vx,T x }. (25) 



We solve this perturbatively in v can , choosing the transformation so that h is independent of A. Also, we 
demand that T\ and v\ do not contain any small energy denominators. Thus we define: 



dh 
dX 
dv\ 
dX 



0, 

[v\,T x ] 



(26) 
(27) 



Eq. (27) is a choice such that T\ and consequently v\ do not allow any small energy denominators. These 
additional constraints determine T\ and v\, which are given by the following equations: 



[h,T x 



-df\ T 



[vy,Ty]dX' 



(28) 
(29) 



where fx + f x = 1 . Eqs. ( |28| ) and ( p9| ) follow from Eqs. (p5f)-(p7|) and the boundary condition H R | A ^ A = 
H B . Now we solve Eqs. ( |2"8| ) and ( p9[ ) for 7\ and TJa- Given Eq. (|l6|), we need to determine v\ : and H R 
is known. The solution to Eqs. (^8|) and (|2^ ) is: 



_ _(1) _(2) 

VA = W A + v \ 



(1) (2) 

T A =T A +T x +■ 



r c (2) (3) 

dv. = dv A +dv, +■■ 



(30) 
(31) 
(32) 



where the superscripts imply the respective order in the canonical interaction, v can , and these quantities 
are given by: 



h,T- 



h,T x 



_< 3 > 



(3) 



dfx 

d\ 
i-A 



dA [u A , , T x , \ + 6v A 



A 



_(=) dfx — (i) (i) 

vx ^-/aK ,r A ] 



(1) (2) 

U A' ) T A' 



dA' 
dfx -j 



(2) (1) 

U A' > T A' 



_(3) Of A -J ( X (1) 



(2) (1) 

«A , T A 



(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
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A general form of these effective interactions is: 

oo „a 

/ dA>y,l£ ) ]+*>? . (39) 



1>\ 



for i = 2, 3, — , with ' = v can . 

H R is renormalized by requiring it to satisfy coupling coherence . A coupling coherent Hamiltonian 
satisfies: 

S(X, A)H B (A, eA,mA,c(eA,m A ))S' t (A, A) = i? B (A, e A , m A , c(e A , m\)) , (40) 

with the additional requirement that all dependent couplings (only one is shown in the argument of the 
Hamiltonians for simplicity) vanish when the independent marginal couplings are taken to zero. Note 
that there are only a finite number of independent couplings. This is a highly non-trivial constraint on 
the theory and to date has only been solved perturbatively. In this paper, Eq. (fio|) is solved to order e 2 , 
which turns out to be fairly simple because e does not run until order e 3 . 

Now we write the solution to second order for . From Eq. (|35|) we obtain: 



/ (AA) (AA) > 

E \-f- + -^- (4i) 

k \ lk jk ) 

where = £ dX ' h , jk ^ . (42) 

(2) 

Sv Ai j will be determined (in §111) by requiring the conditions of coupling coherence to be satisfied. 
These previous equations are valid for an arbitrary similarity function, f\. In this work we will use 
f\ij — #(A — |Ajj|) (a convenient choice for doing analytical calculations). This gives: 

9ikj = {fAik-f\ik)@ikj, (43) 
Q ikj ee 6(\A ik \-\A kj \) . (44) 

To complete this section we write the canonical QED Hamiltonian. We start by dividing H can into 
free and interacting parts: 

H can EE h+ V can , (45) 
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where h is given by Eq. Starting with the QED lagrangian (e > 0, i.e. the charge of the electron is 

-e): 



with 



C QED = + p - m)il> 



(46) 



in a fixed gauge, A + = 0, the constrained degrees of freedom are removed explicitly, producing v can - 
For details of the derivation see §IV.A of Ref. We use the two-component representation chosen by 
Zhang and Harindranath |To| ] . Below we write the resulting Hamiltonian completely. The field operator 
expansions and light-front conventions followed in this paper can be found in Appendix A. The canonical 
Hamiltonian is 



run — H can — h ~\~ V C an i 



cPx^dx 



(47) 
(48) 



where 



'tint — 'Lee'v 



int — i**ee*y i flecyy tH eeee ) 



(49) 



and 



n eei = eg { -2(— • A ) + a • A x — + a ■ A x } £ , 



0+ 



H 



eeyy 



WW 



(50) 

(51) 
(52) 



Note: i — 1, 2 only; e.g., a ■ d x = & L d l = a 1 {—d\) + er 2 (— c?2); {a 1 } are the standard 2 x 2 Pauli matrices. 
Also, /i is given by Eq. (p^). 



B. Step two: diagonalization of H F 



The second step in the similarity hamiltonian approach is to solve for the spectrum of H R . The 
Schrodinger equation for eigenstates of H R is: 

9 



52(i\H*\3){3\**.H<P)) = E N {i\* RiN (P)) . (53) 

3 

See Eq. (|l^) and the comment immediately following it for an explanation of the notation. 'N' labels all 
the quantum numbers of the state, and is discrete for bound states and continuous for scattering states. 
E N = — v+ N , A-P = (2m + B N ) 2 , and "P' is the total momentum of the state of physical interest (for 
this paper, positronium) . 

Solving this eigenvalue equation exactly is not feasible, because all sectors are still coupled: 

\*«An = Ei i )< i i*-«^)) = Ei e ^')>(^')i* i? . JV (p)) + 

i i' 

+ £ |ee 7 (i')><ee 7 (0l* R ,^)> + ]T \eeee(i>))(eeee(i')\* R . N {V)) + ■■■ . (54) 
%' %' 

We divide H R into two pieces: 

h r = n + (H R -n )=n + v , (55) 

diagonalize H exactly, and calculate corrections to the spectrum of H in BSPT with V. The H we 
choose for positronium is: 

H = h + Y,\eni')){ee(i')\V G \e^j'))(ee(j')\. (56) 

i'j' 



h is the free Hamiltonian given in Eq. (|12|). We are assuming photons couple perturbatively to the |ee) 
sector, which must be justified a posteriori. V c is the Coulomb interaction and will be written explicitly 
later. Note that the lowest order low-lying spectrum of the complete H R is identical to that of h + V c as 
long as the limit in Eq. ([!]) is taken. 

We close this section by writing the standard BSPT Raleigh- Schrodinger formulae. For simplicity, we 
write the formulae for the non-degenerate case j|Xl[j : 

(H + V)\* RiN (V)) = E N \* R:N (V)), (57) 
U \^ N {V)) = E N \^ N {V)) , (58) 

N M 



{4> N (V)\v\4, M (V)) 2 



E » = £ » + u (v\ui, <v\\ + 2- ^—f. + ° {v > ' (60) 
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where V is the total three-momentum of the state and "N" labels the total mass of the state. These 
general formulae will be used below in §111 to solve for positronium's spin structure. Note that for the 

V ±2 +M 2 V ±2 +M 2 

light-front case: E N = — v+ N and £ N = ^ — — . 



III. POSITRONIUM'S SPIN STRUCTURE 

Now we will apply the procedure outlined in §11 to obtain positronium's ground state spin splittings 
to order a 4 . We will now give a brief overview of this section. First, we derive H R to second order 
in e. This includes a discussion of the effective fermion self-energy, but the photon self-energy and 
electromagnetic coupling do not run at this order. Then we move on to the diagonalization of H R . This 
starts with a discussion of our zeroth order Hamiltonian, 7i , which will be treated nonperturbatively. 
This includes a discussion of a coordinate change that takes (x S [0,1]) — > (n z £ [— oo,oo]), which 
allows easier identification of H - We solve for the spectrum of TL exactly, which among other things, 
fixes the a-scaling of the momenta in the matrix elements in BSPT. Then we move on to a derivation of 
the perturbative effects coming from low-energy (energy transfer below A) photon emission, absorption 
and annihilation at order e 2 , which includes a discussion of the full electron and positron self-energies 
and a derivation to order e 2 of the complete exchange and annihilation interactions. Given this, we 
determine the range of A that allows the effects of low-energy (energy transfer below A) photon emission 
and absorption to be transferred to the effective interactions in the |ee) sector alone, and at the same 
time, docs not cut into the nonperturbative features of the solutions of TL . Finally, we proceed with 
BSPT in V noting that all shifts appear in the few-body sector, |ee), alone. 

A. Derivation of H R to second order 

From §11, the final renormalized Hamiltonian to second order is given by: 

r / (AA) (AA)\ ^ 

m R \j) = hii K + (Vcanh, + E (V can ) ik (v can ) kj ( ^gL + E^- \ + + {^) \ . (61) 

9ikj is § iven in Ec L- © and "™» is S iven in Ec L- ©■ 
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1. Renormalization issues 



The form of Sv A follows from the constraint that H R satisfies coupling coherence. To order e 2 the 

fermion and photon masses run, but the coupling does not. First, we discuss the result for the electron 

self-energy coming from the second-order effective interactions in H R . We skip the tedious but simple 

details of the calculation of the matrix element, but the interested reader should consult Appendix A, 

where we have collected our light-front conventions for our field expansions and commutation relations. 

We use the rules in Ref. Jl3| to calculate our matrix elements in this work. Specifically, we calculate a free 

matrix element of H R given in Eq. ( |3l|) in the electron self-energy channel. The results are summarized 

by the following equation which also shows how the running electron mass squared, m 2 , is defined]^] 

K 2 +m\ (e(3)e(4)|gJe(l)e(2))| 3ei/ _ rTO P+ 2 K * + m * 

x (e(3)e(4)|e(l)e(2)) l-x 

- ^ + -<<'-^'> + fe *' +C V), (62) 

X X 



where 



(2) a ( 3X' 2 x + m 2 rn (— + \ x) m 2 + \' 2 x 

\> = — — I H 5 - 3m 2 log 5 

A 2tt 1 2 m 2 + \' 2 X & \ m 2 



a 



In the respective ST,^, terms of Eq. (^), A' = A and A. In these formulae, 

Patron =Pi = {xV + ,k + x V^) and V poaitronium = V = {V + ,V ± ) , 

and |e(l)) ( or |e(l)) ) is a state of the free Hamiltonian, h, with spin and momenta coordinates labeled 
by "1." Note that we are forced to introduce an infrared regulator, e. This is introduced by requiring all 
longitudinal momenta (real, virtual or instantaneous) to satisfy: 

l ^ ^ e = +< ( 64 ) 

where V + is the total longitudinal momentum of the physical state. The absolute value sign is required 
for instantaneous lines. Physical results can not depend on this infrared regulator, e, and in this QED 

6 Note that the energy dependence on the electron's relative transverse momentum, re, does not change. 
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calculation we show that treating the photon perturbatively leads to an exact cancellation of this infrared 
divergence in the above running electron mass squared, m\. 

We constrain the electron mass to run coherently with the cutoff which from Eq. (|io| ) and Eq. ( |62|) 
amounts to the requirement: 



This fixes the mass counterterm: 7 



Sv A 



O e 4 



+ Sv ( l ) +0(e i ) 



A^X 



(65) 



Sv? = S^+0(e 4 ) 



(66) 



and to second order the fermion mass renormalization is complete. 

For arbitrary A, the photon mass also runs at order e 2 . The discussion follows that of the electron 
mass except for the fact that the running photon mass is infrared finite. For A 2 < (2m) 2 , the photon 
mass does not run because pair production is no longer possible. Thus, for A 2 <C m 2 the photon mass is 
zero to all orders in perturbation theory. There are additional difficulties with marginal operators that 
are encountered at C(e 3 ), but this is beyond the focus of this paper. 

2. H R to order e 2 : exchange and annihilation channels 



To complete the derivation of H R to second order we need to write the coherent interactions for 
the exchange and annihilation channels in the |ee) sector. At second order, these come from tree level 
diagrams, with no divergences or running couplings, thus the results follow from: 



Sv 



(2) 



[v x , , T x , \d\ 



To show that 5v A produces a coherent interaction recall Eq. We have: 



/ [v$,7$]dX- [,l?,T^]dA' + 0(e 3 ) 

J A JX 



,<.,,„ - J [v x ,\T y '}dX' + 0(e 3 ) 



(67) 



(68) 



7 Actually, any finite C(e 2 ) scale independent term could be added to the counterterm, and Eq. (^) would still be 
satisfied. However, higher order renormalization reveals that this term would not be coherent, and so is excluded. 
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(69) 



A— »A 



which satisfies the coupling coherence constraint, Eq. ([Io|). At second order this seems trivial, but at 
higher orders the constraint that only e and m run independently with the cutoff places severe constraints 
on the Hamiltonian. 

Given this second order interaction, the free matrix elements of H R , shown in Eq. ( |6l| ) , in the exchange 
and annihilation channels are: 
Exchange Channel 

<e(3)e(4)|#Je(l)e(2))L c _ 



X, exchange 



16tt 3 8 3 (pi +P2-P3 - Pi)\/xx'(l - x)(l — x') 
Vi + V 2 +0{e A ) , 



(70) 



where 



Vi = - e 2 iVi 6 (A 2 



Ml - M' a 



^(|A 1 |-|A 2 |)^(|A 1 |-A 2 ) e(|A 2 |-|A 1 |)^(|A 2 |-A 2 



DEN! 



DEN 2 



V 2 = 



(A 2 -Ml- M' ' 



(x - x>)2 S1S3 S2Si 
The variables are defined as follows (see Figure 1 also): 

• Pi = (xV + ,k + xV^) , p 2 = ((1 - x)V + ,-k+ (1 - x)V^) 

• P3 = (x'V + , k' + x'V 1 -) , p A = ((1 - x')V + , -k' + (1 - x')V~ 

• Ni = 6 slS3 5 S2Si Ti ■ T% - 2m 2 (5 S2 s 4 <5 S2 s 1 <S, 

+ imV2(x' - x) ( -^l&aiaa^sta^ ' T{ 



Si 
XX 1 

A „.i' l \ 



S S1 xx'(l- x)(l- x') 

S2 



T? = - 



• Ai 

• DENi 
• X 2 



(l-a;)(l-x')^ 4S2,5siS3 ^ 2 

x — x' 1 — x' 1 — x ' 1 x — x' x' x 

n l +i s tij K j (s = ± 1 and s = — s) ; ei 2 = — e 2 i = 1 , en = e 22 = 
DENi DEN 2 

— ' A 2 = 

x' — x x' — X 

(/tar' — k'x) 2 (mi — mx'Y 



XX' 



XX' 



, DEN 2 = DEN 1 \ x ^ XtX ,^ 



k 2 + m 2 , ,, 2 k + m 
A4 = — 



x(l — x) 



x'(l-x') 



(71) 
(72) 
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Annihilation Channel 



V, 



(e(3)e(4)|if R |e(l)e(2))L 



X^annihil 



16tt 3 5 3 (pi -\- P2 — P3 — P4)\/xx'(1 — x)(l — x f ) 
= V 3 + V 4 + (D(e i ), 



(73) 



where 



and 



V 3 = e 2 N 2 9 (A 2 



Ml-M 



,2 



'e(M 2 -M' 2 )e(M 2 -\ 2 ) 



Mi 



+ 



e (m' q 2 Ml) e (m' q 2 a 2 ) 



M 



/2 



V 4 = 4e 2 6 (A 2 - \M 2 - M' Q 2 \) S Bll3 5 a3lt , 



2m 2 



xx'(l - x)(l - x 1 ) 



1V2 



si 



x(l — x) 



S 4 



x'(l - X 1 ) 



Ti = 



K' 1 (S 3 ) k'\s 3 ) rjri _ «*(«!) K 4 (si) 



1-a;' 



1-a; 



(74) 
(75) 



V 2 and V4 are canonical instantaneous exchange and annihilation interactions, respectively, with widths 
restricted by the regulating function, f\. V\ and V 3 are effective interactions that arise because photon 
emission and annihilation have vertices with widths restricted by the regulating function, f\. 

B. Diagonalization of H R 

First we discuss the lowest order spectrum of H R , after which we discuss BSPT, renormalization 
and a limiting procedure which allows the effects of low-energy (energy transfer below A) emission to be 
transferred to the |ee) sector alone. 

1. 7i , a coordinate change and its exact spectrum 



H in the |ee) sector is 



15 



Ho = h + V c , (76) 



where h is the free Hamiltonian given in Eq. (|12|), and V c is given by (using the same variables defined 

below Eq. (f7|); note, k z is defined below by Eq. @) 

v = (e(3)e(4)|V c |e(l)e(2)) ^_ 16m 2 e 2 ^ lS3 ^ 4 

16^ 3 5 3 (pi +p 2 -pa - Pi )^xx'{l - x)(l - x') (k-k') 2 + {k z -k' z ) 2 ' 

In all other sectors we choose TL = h. TL Q in the |ee) sector was motivated from the form of our second 

order renormalized Hamiltonian, H R , and arises from a nonrelativistic limit of the instantaneous photon 

exchange interaction combined with the two time orderings of the dynamical photon exchange interaction. 

We choose it to simplify positronium bound-state calculations. Other choices are possible, and must be 

used to study problems such as photon emission. Later, in BSPT this choice is shown to produce the 

leading order contribution to positronium's mass as long as the limit in Eq ([!]) is taken. 

The coordinates k z and k' z in Eq. ( [77j ) follow from a standard coordinate transformation that takes 

the range of longitudinal momentum fraction, x S [0, 1] to k z S [— oo, oo]. This coordinate change is: 

' ' (78) 



2 2y/K 2 + K 2 z +m 2 
We introduce a new three-vector defined as: 

P = (/S, Kz) ■ (79) 

Note that 

Ml ^ ^f±^=4(m 2 + p 2 ) (80) 
a;(l — x) 

is invariant with respect to rotations in the space of vectors p. The nonrelativistic assumption of Eqs. (||) 
and (|jj|) in terms of this three-vector becomes: 

IpI 



ni 



0{a) . (81) 



Note the simple forms that our "exchange channel denominators" take in the nonrelativistic limit 

px5 



/n9 (k, — k')(x> 2 — p' 2 ) 
DEN l = p - p') 2 - — — — + O 

m 

DEN 2 = (p-p0 2 + (^-<)(P 2 -P' 2 ) +0 



1G 



TO 2 



P 



5 



m 2 



(82) 
(83) 



Also note the form that the longitudinal momentum fraction transferred between the electron and positron 
takes: 



x — x 



2m 



4m 3 



O 



(84) 



These formulae are used throughout this paper. 

Now we describe the leading Schrodinger equation. We seek solutions of the following eigenvalue 
equation: 



H |Vjv(P)) = £n\i}>n(V)) 



(85) 



where £n = — ~- is diagonal with respect to the different particle sectors, thus we can solve 



Eq. (\8q) sector by sector. In all sectors other than |ee), Tt = h, and the solution is trivial. For the |ee) 
sector, a general \iPn(P)) is: 

\i>N{V)) = E / \fptvt^5\V - pi - P2 ) 4>n(xks iS2 ) bl (pO dl ( P2 ) |0) , (86) 



with norm: 



(4>n(PMn>(P')) = 6 NN ,16ir 3 V+6 3 (V - V) 

EJ d 2 K Ldx ~ 
-— - (f) N (xKSxS2)(f>N'{XKSiS2) = 6 NN ' 
WTT A 



The tilde on 4>n will be notationally convenient below. In the |ee) sector, Eq. (B5J) becomes: 



Ml 



/2 , 2 

k + m 



x'(l - X 1 ) 

After the above coordinate change, this becomes: 



J ^ 2 K In 

[X K S3S4) = 2^ ^3 V c (f) N (XKSiS2) ■ 



M 2 N ~ 4(m 2 + p' 2 )) (Mp'^) = E / 



d 3 p^J(p)J(p>) 



16tt 3 



V c 0Af(pSiS 2 ) 



(87) 



where the tilde on the wavefunction has been removed by redefining the norm in a convenient fashion: 

0~NN' = ^ - 16 ^° (/)* N (xKSlS 2 )(f>N>(xKSlS2) = ^ / d 3 p^^(f)* N (pSlS2)(/>N'{pSlS2) 

" ' (89) 



^2 J d3 P 0Af(pSlS2)0Af (PS1S2) 
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and the Jacobian of the transformation is given by: 



. . dx k 2 + m 2 . „. 

Jp = -t- = — rrr. 90 

«K 2 2(p 2 +TO 2 )2 



Note that the Jacobian factor in Eq. (|38|) satisfies: 

«^('- P ' + 2 ^f + ^ ^(^^---))- 1911 

Before defining 7i in the |ee) sector we mention a subtle but important point in the definition of Ti . 
TL in the \ee) sector will not be defined by Eq. §8^). Rather, it will be defined by taking the leading 
order nonrelativistic expansion of the Jacobian factor in Eq. ( jS8|) . This gives 

(M 2 N - 4(m 2 + p' 2 )) Mp's 3 s 4 ) = E / ^itffi V c <Pn(psi8 2 ) , (92) 

where V c is defined in Eq. ( |77| ) . This H will be diagonalized exactly, and the subsequent BSPT will be 
set up as an expansion in V = H R — TL Q . First, we discuss the exact diagonalization of TL D . 
Putting the expression for V c into Eq. (^||) results in the following equation: 



I p' 2 \ a f d 3 p 
-B N H <Mp's 3 s 4 ) = 7T2 I 1 7T2 ( / , W(PS3S 4 ) ■ (93) 

\ 77i y 1-k 1 j (p- p') 

This is recognized as the familiar nonrelativistic Schrodinger equation for positronium. Note that we 
have defined a leading order binding energy Bn, as: 

M 2 N = Am 2 + AmB N . (94) 

Note the difference in the definition of this leading order binding energy and the full binding energy as 
given by M 2 = (2m + B N ) 2 (see Appendix D for further discussion of this difference). 

To proceed with the solution of Eq. ( |93| ) note that there is no spin dependence in the operator so the 
spin part just factors out: 

<Pfi,s B ,s_(p's3S4) = <Pfi(p')8 Se s 3 5s_S4 ■ (95) 

We rewrote N as (p, s e , s_), where (s e , s ) label the spin quantum numbers and [i labels all other quantum 
numbers, which are discrete for the bound states and continuous for the scattering states. 
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The solutions to Eq. (£)3j) are well known. For B N < 0, following Fock [[12| , we change coordinates 
according to 



(uo,u) 
cos(w) 

p . 



4-p 2 
4 + p 2 



sin(oj) = sin (a;) (sin(#) cos(0), sin(9) sin(</>), cos(#)) 
2e„p 



el + p 2 



(96) 
(97) 
(98) 

(99) 



Useful relations implied by this coordinate change are in Appendix B. Note that in our notation we 
anticipate that \i will be given by (n, Z, m l ), the usual principal and angular momentum quantum numbers, 
and that the leading order binding will depend only on the principal quantum number, n. Given this, 
Eq. ( p3| ) becomes 



MM) = 



2tt 2 2e n J \u- u'\ 



(100) 



where 



-5— <Mp) • 



4(0- 



(101) 



Using Eq. (182) of Appendix B, this is seen to have the following solution: 



,„s am 2ir 2 

</; M (fi) = Y^fi) with ___ = ! 



(102) 



where Y^(O) is a hyperspherical harmonic. Thus, 



ma ma 
en = _ and B N = -^ 



(103) 



This is the standard nonrelativistic solution for the bound states of positronium to order a 2 . This 
completes the solution of 7i for the bound states. The scattering \ee) states are also needed in our 
BSPT calculation. We use propagator techniques to include these scattering states where required (see 
Appendix C). 
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2. BSPT, renormalization and a limit 



Here we use the BSPT formulae (appropriately generalized to the degenerate case) of §11. B to analyze 
positronium's spin structure. The potential to be used in BSPT is: 



V 



(104) 



where the eigenvalue equation for H is given by Eq. (J92|), and H R to second order is given in §111. A. We 
will be perturbing about the nonpcrturbative eigenstates of TL . 

First, we discuss electron mass renormalization. In second order BSPT there is an electron mass shift 
coming from the f\v carl part of H R , with v can given by J d 2 x A ~dx~Ti. ee -y (see Eq. fl5C|)). This is photon 
emission and absorption restricted by the regulating function, f\. The calculation is similar to that of 
§111. A. 1. Assuming (M 2 — M. 2 ) — 0(e 2 ), this electron mass-squared shift is 



5m 



(105) 



<5X A is the same function that was defined in Eq. (|j§. Using this result (Eq. (gog)) one obtains: 
(M. 2 , — M. 2 ) = (4m 2 + AmB N — 4(m 2 + p 2 )) = 0(e 4 ), and our initial assumption is satisfied. When this 
is combined with the only other second order electron mass shift, to 2 , of Eq. (|6^) we have, for the full 
electron mass-squared, to second order, 



to 2 + dm 2 

m 2 - (W 2) 
m 2 + C(e 4 ) . 



bv K 



-6K 



0(e 4 ) 



(106) 



In this last step we recalled the result from Eq. (|6g). We see that to second order, the full electron mass 
is given by the electron mass in the free Hamiltonian, h. Also, as promised below Eq. (Q), we see that 
treating photons perturbatively has led to an exact cancellation of the infrared divergence in the running 
mass, m 2 : the full electron mass, to 2 , to second order is infrared finite. 

Now we move on to the discussion of BSPT. The only channels to order e 2 are exchange and annihila- 
tion. Parts of these effective interactions are given in §111. A. 2 . We also need to include the perturbative 
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mixing of the \eerf) and I7) sectors with the |ee) sector arising from f\v can , with v can — J d 2 x ± dx~Tt eei . 
In second order BSPT this gives rise to the following effective interactions that must be added to V x elcchange 
and V x annihil of Eqs. ( f70| ) and (|73| ) respectively. 
Exchange Channel 

- e ^ 1 fl(A 2 -|A 1 |)g(A^-|A a |) 
Vb DEN, ' (1 ° 7) 

with 

DEN 3 = ( K - K ') 2 + i(x-x')A+\x-x'\ (^(m 2 +M i2 ) - M] 

and 

k 2 + m 2 k' 2 + m 2 k' 2 + m 2 k 2 + rn 2 

A = 1 1 ~ + ; • 

1 — x 1 — x x x 



Annihilation Channel 



(A 2 - M 2 Q ) e (a 2 - M' a 



V 6 = e 2 N 2 . (108) 

Note that in a nonrelativistic expansion (after the coordinate change of Eq. (|7|)), the above "exchange 
channel denominator" becomes 



DEN 3 = (p-p>) 2 + \ x - x >\(^-(m 2 +M> 2 )-M 2 n ^ + 
The full exchange and annihilation channel interactions to order e 2 are 



O 



P 



— m 



(109) 



V^ change = V^ xehange + V 5 , (110) 
V annihil = V X:annihU +V 6 , (111) 

where Eqs. ( f70[) and (|7^ ) give V x cxcha „ 3C and V x annihil respectively. 

One way to summarize the results, recalling the form of Eq. (^8|) an d the norm in Eq. (jsi^), is to state: 
the full order e 2 effective interactions give rise to the following first order BSPT shift of the bound-state 
mass-squared spectrum of 7i : 
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$ M 2 (s 3 ,S 4 ;s 1 ,S 2 ) = {<f>n,l,m,i,a a ,Si\V\<f> n ,l tmi , sl ,a a ) 

d 3 p d 3 p'<j)* , (p')V(p', s 3 , s 4 ; p, si, s 2 )<£„,2, mi (p) , (112) 



where 



The Dirac notation in Eq. (112) will be used in the remainder of this paper. See Eqs. (f77|), (110) and ( [111 



for V c , V exchan e and V annihit respectively. The interaction V must be diagonalized in the degenerate spin 

space following the standard rules of degenerate BSPT. Note that V needs to be considered in second 

order BSPT in this paper also. 

The diagonalization of V in the degenerate spin space follows shortly, but first recall from §11 the 

range of A that allows the effects of low-energy (energy transfer below A) photon emission and absorption 

to be transferred to the effective interactions in the |ee) sector alone, and at the same time does not 

remove the nonperturbative bound-state physics of interest: 

|M 2 -(2m) 2 | A 2 , N 

' V+ ' « « W • ( 114 ) 

After the solutions of TL Q are known the a-scaling in all BSPT matrix elements is known and the bounds 



in Eq. (114) become 

m 2 a 2 < A 2 < m 2 a . (115) 



This is satisfied under the following limit: 



to 2 a 2 



A 2 



A 2 — > a fixed number , (H6) 
-> 0, (117) 
-> oo . (118) 



m 2 a 



A 2 



Given the nonrelativistic limit: 



m 
A 2 



a — >• , (119) 

2 



oo , (120) 
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this implies 

^-cxa~, (121) 

where 

2 < k < 4 . (122) 
Note that this "window of opportunity" is available to us because one, we have introduced an adjustable 

\2 , 

similarity hamiltonian energy scale, ^r, into the theory, and two, QED is a theory with two dynamical 

2 2 2 

energy scales, "L° and a fact known for a long time, and the reason that QED calculations have 
been so successful over the years. 

Given the above limit (Eqs. (|Ill|-(|il~8|)), 



A 2 - 4|p 2 - p' 2 |J , 9 (IAxI - A 2 ) , 9 (| A 2 | - A 2 ) — 1 (123) 

(4(p 2 + m 2 ) - A 2 ) , 9 (a( P ' 2 + m 2 ) - A 2 ) — > 1 (124) 

0(A 2 -|A 1 |),0(A 2 -|A 2 |) — >0 (125) 

9 (A 2 - 4(p 2 + to 2 )) , (A 2 - 4(p' 2 + to 2 )) — > . (126) 



Now we proceed with the diagonalization of V in the degenerate spin space (see Eqs. (112) and (113)). 
We will calculate all corrections to order a 4 that arise in the spin splitting structure of the ground state 



spectrum of TL via BSPT. First, we write V more explicitly given the above limits in Eqs. (123)-(126): 



107T 3 Zm \ 4to z \ mr 



e 2 N, 4e 2 r r , e 2 N 2 2 

+ 4e d Sl s 2 d s 



D^iV 4 (x - x 1 ) 2 Us ^ Us ^ T DEN 5 
11 , 



167r 3 2m 



where 



DEN4 DENi DEN 2 

! _0(X 2 -M^ 2 ) 9(m' 2 -M 2 
DEN 5 = 7W 2 + 7^2 
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12 = (9 (DENi - DEN 2 ) (128) 

(129) 



Note that we have expanded out the Jacobian factors as given by Eq. (J9lj). Also, DEN\ and DEN2 are 
defined below Eq. (|7^) and written in their expanded version in Eqs. (^) and (^3|) respectively. Finally, 
iVi and N2 are written below Eqs. ((72J) and (|75| ) respectively. 

Now, since the eigenstate wavefunctions of 7i force p to scale as p ~ ma, it is useful to note the 
a-scaling of specific terms in V. Recalling that we are always assuming a — ► (without which our 
matrix elements would not have a well-defined a scaling), we see the following structure arising 

V = v w +V W +V (2) +-.. , (130) 

where V ' scales as V ( ~ a s . Thus in first-order BSPT these respective terms contribute 

5 m M 2 NN , = (0 N \V lS) \cj) N ,) ~ a 3+s . (131) 

So, to be consistent to order a 4 we need to look at all the matrix elements V ' with S < l. 8 
Before proceeding to write out these expressions for V , we note the following facts: 

• In this work we will only calculate the spin splittings, so any constants along the diagonal in spin 
space do not contribute. 

• Since we are only working to order a , obviously any splittings that can be shown to contribute at 
order a 4+k with k > need not be calculated. 

• Symmetries of the integrand can be used to simplify expressions immensely. 

One final discussion that we must have, before we write out these expressions for V , is how we are 
going to deal with DEN4 and DEN$ defined above.[] These denominators are dealt with by noting the 
following formulae: 

9{a-b) 0(6 - a) 1 6(a -b)+ 0(b - a) 1 6 (a — b) + 6{b — a) 

a + 6 ~ 2 a + 2 6 

1 6(a - b) - 6(b - a) 1 6 (a — b) — 6{b — a) 
+ 2 a 2 6 



For example, . e p . ,.> ~ => S = . 
(p-p ) « 



9 Actually the DEN5 term is handled with analogous techniques as the DEN4 term, and has even smaller corrections 
than those of DEN4. Thus, we will just discuss the DEN4 term in what follows and here state the result for the DEN5 
term: Take DEN$ > 4m 2 ; the corrections to this start shifting the bound state mass at order a 6 . 
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+ l -{6(a-b)-6(b-a)) (I -I 



l\a-b\ (I 
2 a — b \a 
l\a-b\ 
2 a b ' 



(132) 



To proceed it is useful to note 

DEN l 

DEN 1 
1 1 



(p-p') - 

DEN 2 - 
1 



(^ 2 -4)(p 2 -p' 2 ) 

m 

2( Kz -4)(p 2 -p' 2 ) , 



o 



o 



m 



( £ ) 5 
(9 



771 



p\ 2 1 



(133) 
(134) 
(135) 



2 \DENi ' DEN 2 J (p - p') 2 m 2 (p-p') 6 

Especially note that this last equation scales as: + 1 4- a 2 + ■ ■ •, i.e. the corrections start at order 1 
(not order i); this implies that only the ( p _p/)ii term of Eq. (135) contributes to the spin splittings to 



order or. But we still have to discuss the second term that arises in Eq. (132). This term is given by 



DEN 4 



1 \DENx - DEN 2 



2(K z -«;)(p 2 -p' 2 ) 



O(a ) 



(136) 



2 DEN 1 DEN 2 2 ( p _ p ') 4 

Including N%, this starts out as an 0{a?) spin conserving contribution. The next order contribution is 
C(a 4 ) with spin structure, but is odd under p <-> p', and thus integrates to zero. However, the 0(a 3 ) 
spin conserving term appears to lead to an order a 4 shift to the spin splittings in second order BSPT 
when the cross terms with V (0) of Eq. (138) are considered;]^] however, these cross term contributions 
add to zero due to the facts that the 0{a 3 ) term including Eq. (U3q) conserves spin and the 0(a 3 ) term 



including Eq. (13£) is even while the term from Eq. (138) is odd under p <-> p'. 

To summarize the preceding discussion of DEN4 and DEN5, we can say that through order a 4 , for 
the spin splittings of positronium, there are no relativistic corrections to the following replacements: 



DENi — > (p - p') 2 and DEN 5 



4m 2 



(137) 



This is valid for the ground and excited states, but in what follows we specialize to the ground state for 
simplicity. 



10t/(°' 



V of Eq. (h_3§) comes from the first term of Eq. (|135|) combined with the complete next to leading order term in JVi . 
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Given this general conclusion about DEN4 and DEN 5 , we list the pieces of V that contribute to 
positronium's ground state spin splittings through order a 4 . Explicitly, as far as the a-scaling goes, we 
need to consider 

V (p's 3 s 4 ;psis 2 ) = - — - I w (0) (p's 3 S4;psis 2 ) , (138) 

where 

w (0) (p's 3 s4;psis 2 ) = {5 Sl s 3 5 S2S4 fi(p's 3 S4;psiS2) + 6 SlS3 5 S2 s 4 f 2 (p's 3 S4;psiS2)) , (139) 
/i(p's 3 s4;psis 2 ) = s^Ky - k' v ) - i(K x - k' x ) , (140) 
/2(p's 3 S4;psis 2 ) = S4,(k v - n' y ) + i{n x - n' x ) . (141) 

Recall that Si = ±1 (i = 1, 2, 3, 4) only. The only other interaction we need to consider is 

2 

V (p's3S4;psis 2 ) = (c an d SlS2 S SlS J S3Si + c ex S S2l J S2ll S S3ll + 



1 _ (k-k 1 ) 2 

2 C " (p - p' 



+ ( - c„ L . 2 5 S1?2 <5 S3?4 ) . (142) 



The constants c ex and c an were introduced only to distinguish the terms that arise from the 'exchange' 
and 'annihilation' channels respectively, and c ex = c an = 1 will be used in the remainder of this work. 

Two simplifications were made in deriving V W . First, we did not include terms that are a constant 
along the diagonal in spin space, because these do not contribute to the spin splittings to order a 4 . 
Second, we noted that terms of the following type integrate to zero: 

C ihirpKj K/vKi'y. /? X K ) 

(4W 3 .J (p_p^2 -Ko,o, sl ,J = (0,0,0) . (143) 

and thus were not included in the definition of V { ' . 

The ground state spin splitting to order a 4 contains contributions from V a> in first order BSPT 
(V (0) vanishes in first order BSPT) and v'' 0> in second order BSPT. We begin with the first order BSPT 
calculation. These results are shown in Figure 2. Then we perform the second order BSPT calculation. 
The combined results of first and second order BSPT are shown in Figure 3. 
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First Order BSPT : 

The lowest order wavefunctions are given near the end of §111. B.l (see Appendix B for the hyper- 
spherical harmonics). V W in first order BSPT contributes the following to positronium's ground state 
mass squared: 

8Ml = S il) M 2 {s 3 , S4 ; Sl , S2 ) 

/l 1 1 

d "P rf V (e 2 +p2)2 (e 2 +p/ 2 )2 ^ (1) (P'- «3, S4-, P, Sl,S 2 ) , (144) 



where 



8e? , ma 
N = — - and ei = — — 

7H 2n 



n=l 



(145) 



Using the rotational symmetry of the integrand we can make the substitution 



(146) 



(p-p') 2 (p-p') 2 3' 

After this, the remaining integrals are trivial and the splittings that arise from diagonalization of the 
8M 2 matrix in spin space are: 

(l\8M 2 \l) = -m 2 a 4 , (147) 

(2\5M*\2) = \m 2 a 4 , (148) 

(3\6M 2 \3) = m 2 a 4 , (149) 

(A\SM 2 \A) = m 2 a 4 , (150) 

where 

{|l) S ' + ->^- + > ,| 2 ) S ' + -» + '- + » , l 3) s| -),|4) S | + +) }. (151) 

Figure 2 shows these results, which taken alone do not produce the degeneracies required by rotational 

invariance. 

Second Order BSPT : 

V {0> gives rise to the following contribution to positronium's ground state mass squared in second 
order BSPT: 
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SM 2 = S^M 2 (s 3 ,s 4 ; Sl ,s 2 ) 

,0,0,si,s 3 / (1KO\ 

- 2^ P M\- Ml • (152) 

s c: s FM ^(l,o,0) 1 n 

Recall that /i = (n,l,mi), the usual principal and angular momentum quantum numbers of nonrelativis- 
tic positronium. The calculation of SM 2 is tedious, but can be done analytically. This calculation is 



performed in Appendix C. The result is (see Eq. (248) in Appendix C): 



2 4 

777 / V 

5Ml = -^ r (3g 1 +g 2 ), (153) 



where <?i and c/2 are given in Appendix C in Eqs. ( J21C| ) and ( pi 1[ ) respectively. 

Now we combine the 5M 2 and SM 2 matrices and diagonalize the result. The combined matrix is 
given by: 

SM 2 + SM 2 



-^(3.gi+52). (154) 



The eigenvalues are: 



(l\SM 2 + SM 2 \l) = -^m 2 a\ (155) 



(2\SM 2 + SM 2 \2) = \m 2 a i , (156) 
(3\SM 2 + SM 2 \3) = |m 2 a 4 , (157) 
<4|(5A/ 1 2 + <5M 2 2 | 4) = \m 2 a i , (158) 



and the corresponding eigenvectors are: 

> 5 ' + ->^- + > ,12) 5 l + -> + |- + » ,|3) 5 |-), l 4) 5 | + + )}. (1 59) 

Figure 3 displays these results. These results lead to the well-known result, ^a 2 Ryd, as detailed in 
Appendix D. We see rotational invariance in the degeneracies of the ground state n = 1 levels exactly 
maintained at order a 4 . 
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IV. REMAINING CORRECTIONS TO THE POSITRONIUM SPECTRUM 



We have shown how one obtains spin splittings at 0(a 4 ), starting with a light-front Hamiltonian 
renormalized through 0(a). Even though the spin splittings are correct through 0(a 4 ), there are cor- 
rections starting at 0(a 3 ) in the actual position of the levels computed using H R renormalized to 0(a) 
and m 2 a 2 <C A 2 <C m 2 a. We will outline these corrections and argue that higher order counterterms, 
obtained from renormalizing the Hamiltonian through 0(a 2 ) and beyond, fix any errors appearing in the 
physical spectrum computed with H R renormalized to 0(a). 

The dressing of electrons by two photon loops in all allowed time-orderings shift the electron self- 
energy by 0(a 2 \ 2 log 2 (l/e)). In an analytic calculation this contribution can be cancelled by a one-body 
counterterm that is exactly what one expects from higher order renormalization when the Hamiltonian 
is computed to 0(a 2 ). The infrared divergences must cancel leaving a free electron energy p . 

The leading correction to the spectrum we have computed in this paper may be 0(a 2 A 2 log 2 (l/e)), 
coming from the above 0(a 2 ) correction to the electron self-energy. However, with H R renormalized 
through 0(a), there are indications that all the infrared divergences arising in BSPT cancel due to the 
fact that large wavelength photons decouple from neutral positronium states , leading only to residual 
long-range two-body potentials that may still correct the spectrum. Assuming this holds at higher orders 
too, the leading correction from higher order self-energy shifts is 0(a 2 A 2 ). With the cutoffs used above 
the corrections to the electron self-energy can depend on the longitudinal momentum fraction of the 
electron, while if cutoffs involving the momenta of the connected pieces of the matrix elements are used, 
the relevant counterterm is a constant. 



From Eq. (121) in the paper, the constraint placed on A is 

A 2 oc m 2 a* , (160) 

where 2 < k < 4. So the leading order shift to the spectrum arising from higher order self-energy 
corrections is 

5M 2 X ~ m 2 a 2 a* . (161) 
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Another correction to the spectrum from the 0(a) Hamiltonian arises from the interaction containing 
the leading order term of Ni and the term from Eq. Jl3G| ). This shifts the spectrum an amount 

SM 2 2 ~ m 2 a 3 . (162) 

So, after considering the possible corrections, we see that our choice of 7i has produced eigenvalues 
that are correct through (D(m 2 a 2 ). That is, the nonperturbative physics was correctly chosen and is 
not changing. However, fourth order perturbative effective interactions are required to obtain O(m 2 ot ) 
precision for the actual position of the levels — the spin splittings are correct through order a as detailed 
in the paper. 

To be complete, we must also consider coupling and photon mass renormalization. However, there 
are no further corrections to the photon mass if A < 2m because pair production is no longer possible. 
This also affects coupling renormalization. If A < 2rn, no electron loop appears to modify the electron- 
photon vertex. Since this contribution is entirely responsible for renormalization of the marginal part 
of this operator (i.e., the running coupling), we do not expect the charge to renormalize. There will be 
modifications to the irrelevant part of this vertex, but they should not contribute to leading order. 

There are additional counterterms that are not typically encountered in a lagrangian calculation. For 
example, there are tree level counterterms that modify the electron-photon and positron-photon coupling 
at C(e 3 ). These counterterms should produce cancellations with photon exchange contributions to the 
mass, but they must be present for this to occur. Such counterterms may be necessary for all the infrared 
divergences to cancel. 

We offer no full solution to these problems other than renormalizing the Hamiltonian to higher orders 
using a similarity transformation and coupling coherence, as outlined and illustrated in earlier sections. 

V. SUMMARY AND DISCUSSION 

If H R is derived approximately to a finite order in a, to what order in a is the nonperturbative 
spectrum correct? The precise connection between the approximate renormalized Hamiltonian and the 
physical spectrum is only qualitatively understood at present. The results of this paper make it clear that 

30 



in light-front QED to obtain positronium's spin structure through order a 5 , the renormalized Hamilto- 
nian must be derived either through at least second order in a or else via a nonperturbative similarity 
transformation which uses Coulomb states instead of free states in the perturbative expansion for the 
effective interactions. 

We have calculated the order a 4 spin splittings for the n = 1 levels. This restriction leads to a great 
simplification because constants along the diagonal in spin space do not contribute to the splittings. Some 
of the terms that are not calculated because they cancel in the difference are quite complicated and a 
future analytical calculation including these terms would be quite complex, involving complicated sums 
over 9-J symbols for example. Even though we do not calculate these a 4 "radial" shifts in this paper, 
the fact that our spin splittings are correct and can be computed analytically from first principles gives 
us much hope that the procedure outlined in this work is well-defined and at least for QED 3+1 leads to 
consistent bound state calculations. It becomes very interesting to think about the Lamb shift. 
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APPENDIX A: LIGHT FRONT CONVENTIONS 



In this appendix we present our light-front conventions. In this initial paragraph, let A and B be 
arbitrary 4-vectors. An "ET" label implies an equal-time vector; the absence of a label implies a light- 
front vector. Two vectors by definition are related by 

• A^~ = Afirp i Agrp 

• A = Agrp : 2=1,2 

• A" = (A+,A-,A l ) 

• A ET = (A° ET , A ET , A 1 ). 

Then scalars are required to agree by adjusting the metric tensor so that this is so. This fixes the 
light-front metric tensor: 

A^R — n A^R V — 4 M R ET — a ET 4 M R v — A° R° — A 3 R 3 — A 1 R i 
=> 1 = 2.g+- = 2.g_ + = = -.922 • 

Another relevant scalar is of course 



1 = ^: = ^ = -^ = -^ 



The components of the light-front metric tensor not mentioned are zero. Thus, there are factors of two 
in places like 

d 4 x = ^dx+dx-^x 1 - , A_ = g_+A+ = \{A° + A 3 ) = ^(A - A 3 ) , etc. 

Conventionally, x + is chosen to be the light-front time coordinate. This fixes x~ to be the light-front 
longitudinal space coordinate. Also, from the p ■ x scalar: 

P^ = 9+-P~x + + 9-+P + x~ - P l x l = l;P~ x+ + \p +x ~ ~ V %x% > 
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we see that p~ is fixed to be the light-front energy coordinate, and p + is fixed to be the light-front 
longitudinal momentum coordinate. 

All constituents in the forward light-cone in light-front coordinates have p + > 0. This can be seen 
from the following relation (which also shows why particle lines in hamiltonian diagrams are all on mass 
shell): 

d 3 p 



I Wf 2nS{p2 ™W)/(p) = / 

d 2 p ± dp + 6(p+) 



I 



167T 3 p 
/(P) 



fiP) 



167T 3 p+ 

Especially note this last definition of J p , which is a shorthand used in the paper. 
In momentum space the field operators are expanded as (at x + = 0) 

A\x) = J2 [(&.(q)e- iq - x + h.c.), 
Z(x) = E X. I Vp+(bs(p)e- ip - x + 4(p)e +i ^) , 

s=±l Jv 

with 



x 1 = 



V / 



The fermion helicity can only take on the values ±1/2, however we define /13 = s/2; therefore, "s" can only 
take on the values ±1. Note that s = —s. Examples of the commutation (anti-commutation) relations 
and free Fock states are 

[a x (q),a{,(q')} = 16n 3 q+S 3 (q - q')5 xx , , ( 6 3 (p) = 5 2 (p^6(p+) ) , 
{b s (p),bl,(p')} = {d s (p),dl,(p')} = 16n 3 p+6 3 (p-p')d ss , , 
(pisi\p2s 2 ) = 167r 3 p+(5 3 (pi -p 2 )S Sl s 2 , bisi) =6^(pi)|0) , etc. 

The inverse longitudinal derivative can be defined as follows (we define it by putting the momentum 
representations of the field operators in the Hamiltonian, multiplying the terms out explicitly, and then 
replacing the inverse derivative by appropriate factors of longitudinal momentum [x ±i] but nevertheless 
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the following could be used too): 



^Vj/O ) = \ I dy e{x -y )f(y ) 
d 



d+ = 2d- = 2- 



' dx 

<9_e(ar - y-) = 28{x- - y~) , 
e{x) = 6{x)-9{-x) . 

Notice that this is non-local in the longitudinal direction. 
APPENDIX B: HYPERSPHERICAL HARMONICS/FOCK COORDINATE CHANGE 

In this appendix we will list some useful mathematical relations used in the paper. The conventions 



followed in this paper are given in Ref. 1 14 . These hyperspherical harmonics are given by: 



r„(n) = F n ,;, m (fi) = f n j(uj)Yi !m (e,(j)) , (163) 



where 



< \m\ < I < n - 1 . (164) 

These quantum numbers are the standard "hydrogen" quantum numbers. These 3D spherical harmonics, 
Y ltm (0,(j>), are given by JL5J: 



where x = cos (9. These other functions, f n< i(u)), are given by: 



/»,lM = H)S/ 2n( ^ in, 1 ^ 8111 ^ ^ ^ v Cn-i(cosa;), (166) 
w ir(n + l)\ (ctcosw) 1 

where C„_i(cosijj) are Gegenbauer polynomials. For example 

C (y) = 1, C x {y) = 2y, C 2 (y) = 4y 2 - 1 , C 3 (y) = 8y 3 - Ay , 

C 4 {y) = 16y 4 -12y 2 + l, C 5 (y) = 32y 5 - 32y 3 + 6y . (167) 
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The orthogonality and phase relations of these above functions are: 



Yn,l,m — (-1) +mY n,l,-m ' Y l,m ~ ("l)"'^*-™ ! fn,l ~ {—^)fn,l ' ( 168 ) 

<m w = dn = dn (3) du sin 2 u , (169) 

J dQY*Y^ = 5^> , J duo sin 2 u f* tl f n >,i = <W , / dQ, ( ) Y l * m Yi> im i =5u>5 mm > ■ (170) 



The first few hyperspherical harmonics are: 

1 \/2cosw 

'X, 0,0 = , 7; ) -^2.0,0 = 

V2tt 2 tt 

— ie~ 1 ^ sinw sin 6* — i\/2 sin a; cos ( 

i2,X,-X = , J2, X.O = 



7T 

For further harmonics we refer the interested reader to Appendix 2 of Judd's text |l7) , where this is done 
quite nicely 

For the coordinate change in the Coulomb Schrodinger equation (see Eq. (|9^) and the discussion that 
follows it), for B N < 0, we define 

mB N = -e 2 , (172) 
u = (it ,u) , (173) 
u = 008(^ = 4^4, (174) 

u = — sin (a;) = sinfo;) (sin(#) cos(</>), sinlff) sin (</>), cos(#)) 
V 

2e " P (175) 



el + p 2 



Note that Uq + u = 1. Conversely this coordinate change gives: 



We also have: 



P = tx^u, (176) 
1 + u 

4 + = t^-. (177) 
1 + u 



dn p = sm 2 ojdojdn { p 3 ' = ( 2 2e " 2 ") rf 3 p , (178) 

35 



s3 (p-p r ) = f lyL ^p - = (i+ 3 mo)3 ^p - m , (179) 

0<w<7r,0<6»<7r,0<(/)<27r. (180) 



Finally, a most useful relation is: 



IP-P'I 2 = (e " + P ; )( 2 e " + p/) ln-,f . (181) 

This is useful because we can expand \u — u'\ 2 as follows: 

1 ?ir 2 
T—77H = £— (182) 



This completes the discussion of the hyperspherical harmonic mathematical relations used in this work; 
actually, as shown in the next appendix, the SM 2 calculation requires some further formulae, which are 
given as they arc needed. 

APPENDIX C: THE CALCULATION OF SM 2 

i 

In this appendix we perform the following sum analytically: 

5M > = ^ Mj -Ml ■ (183) 

S e ,Se ^(1,0,0) 1 ™ 

Recall that fj, = (n,l,mi), the usual principal and angular momentum quantum numbers of nonrelativistic 
positronium. Also recall that the spin factored completely out of our lowest order Schrodinger equation, 
so to proceed the following notation is useful: 

|^, Se , s _) = |<?V) ® |s e Se) , (184) 
1 = X] I^X^VI ® |s e s-)(s e s-| 

= d3 P |P><Pl ® l s es-)(s e s-| . (185) 

To proceed, define the following Green's function for arbitrary S: 

= z.^— m- (186) 



36 



The factor -J— will turn out to be useful. This Green's function satisfies the familiar Coulomb Green's 

4m 



function equation: 



S 3 (p-p') = (E-*-)G E (p',p) 
m 



a f ,3^ G e(p",p) 



2tt 2 



(p' - p") 5 



where 



and 



p'\G E \p) = G E (p',p) 

E - 4m 2 



E 



4m 



(187) 



(188) 
(189) 



Hostler and Schwinger independently obtained the solution for this Coulomb Green's function in 1964 |18| . 
We find Schwinger 's form useful; the equation he solves is exactly the above equation with the following 
shifts in notation: 



(Ze 



Schwinger 
TFlSchwinger 

E s chwinqer 



a , 

rn 
~2 

E . 



(190) 
(191) 
(192) 



His result is amended because the sum we need has E — Ai\ and does not include fi — (1,0,0). This 
subtraction of the fj, — (1, 0, 0) term amounts to the term "— -t" in Gm below. The details of how this 



arises can be seen in Eqs. 


(|233|)-(|23^) below. With this amendment, Schwinger's result is: 




G^,2(p,p') = 


Gi + Gn + Gm , 




(193) 


Gj = 


<5 3 (P-P') 

E-T ' 






(194) 


Gn = 


a 1 


1 1 




(195) 


2tt 2 e - T ( 


p-p') 2 E-T' ' 




Gin = 


a 4ef 


\l X d P ( 1 IV 


1 


(196) 


27r 2 E-T 


[Jo p WAv-v') 2 + c{i-pY c)_ 


E-T' ' 


where 












T = 


2 

P rp, P ma 

m m 2 




(197) 




C = 


(e 2 + p 2 )( ei +p' 2 ), 




(198) 




E = 


Ml - 4m 2 e 2 
4m m 
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(199) 



The prime on G' M2 denotes the fact that we have subtracted the /i = (1,0,0) part of Ge as required 
by the sum that we have in SM 2 . Note that this Green's function is symmetric under p <-» p' and also 
(jPxtP'x) (PyjPy)i symmetries that will be used in later simplifications of the integrand of SM 2 . 
SM 2 is now: 



SM 2 



/ d 3 pd 3 kd 3 P 'd 3 k'(<p 1 , 0tQ \k)V {0) (k,s 3 ,s 4 ;p,s e ,Se) 

f Gi + Gji + Gni\ jo) i \n/\i \ 

x I — I V (p ,s e ,Se;k ,si,s 2 )(k |01,o,o) 

5M 2 (I) + SM 2 (II) + SM 2 (III) respectively . 



(200) 
(201) 



Now we rewrite this in terms of hyperspherical harmonics and perform the integrations analytically. The 
variables are defined as: 



[u = (u , u)] <-> [ei, p] , [u' = (u' Q , u')] <-> [ei, p'] 
[v = K,v)] «-> [ ei ,k] , [v' = K,v')] «- [ei,k'] . 



(202) 
(203) 



0^ 

See the previous appendix on hyperspherical harmonics for a summary of the mathematical relations that 



we use. The symbols appearing in Eqs. (202) and ( 203 ) are explained there. Note that we use e% in these 
variable definitions, a choice that is completely general and turns out to be useful because we are taking 
expectation values of n = 1 states in this work. The relations we use are: 

4ef 1 



(k'|^i,o,o) = 



'A + k' 2 ) 2 V2^ 



1 



(i + <) 2 



(ej + k' 2 ) 2 4ef 

. d 3 k'= [ei+ a k 3 > <m k , = - , * rfn fc ,. 

8e? (l + <) 3 



(204) 
(205) 
(206) 



Given these, JM 2 becomes: 



SM 2 



ma 



(E — T) (Gi + Gu + Gin) 



32tt 2 J (1 + u )(l + u' a ) 2 

x sJ2^-jy^p) Y ^(^ P 'Ki^k)Y;,(n k ,) , 



(207) 



mm' 



where 
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S = ^w (0) (k,S3,S4;p,s e ,s-) V ( )(p , ,s e ,s-;k , , Sl ,s 2 ) . (208) 



Recall Eq. ( 139 ) for the definition of v^°\ Using the symmetries of the integrand, the sum over spins 
s e and Se can be performed and a simplification is seen to arise. The spin completely factors out of the 
momenta integrations. In other words, we have: 

S = i(3fli+a 2 )(p-p' + k-k'-2p-k') , (209) 
o 

where 

gi = sis 3 + s 2 s 4 , (210) 
.92 = 1 + sis 2 - s 2 s 3 - sis 4 + s 3 s 4 + S1S2S3S4 ■ (211) 

Recall that Sj = ±1, (i — 1,2,3,4); i.e., the '|' has been factored out of these spins.[^] So, in other 
words, instead of having to do sixteen twelve dimensional integrals because the spin and momenta are 
coupled together, we just have to do one twelve dimensional integral that is independent of spin and then 



diagonalize the result in the 4x4 dimensional spin space with the spin dependence given by Eq. ( 20E ) . 
We define the following integral: 

ma. ,„ „. 

X - g^, (212) 

where 



{E-T)(Gi + G n + G n i) 



(i + u ){i + u> y 

x (p-p'+^-2p-k')J2^Y^n p )Y^(n p ,)Y;(n k )Y;,(n k ,) , (213) 



nn 



and 



2 4 

Ttl f\ 

5M 2 2 = --—(3g 1+ g 2 ) X . (214) 

For the quantities £, \ and SM^, the labels I, II and Ii7 imply the respective terms with G/, Gjj and 
Gin above (see Eq. ( |l93| )). Also, the terms a, b and c above correspond to the respective superscripts 
in what follows. This integration will now be performed analytically. 

11 In order to get these simple forms for g\ and 32 and to see this spin/momentum decoupling it was useful to note the 
following simple relation: <5 3s / = ^s(s + s') (true because s 2 = 1). 
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s3 (p-p) = 772 rbw ^^ ~ n ^ = ^ "p') ■ ( 215 ) 



First the three Gj pieces. The mathematical relations used here are: 

MO <>.• 

(ej + P 2)3 

P 2 = 737 — ^ - M °) ' ( 216 ) 

k • k' — ► 3fc z fc; , (217) 

pk' — 3p z fc£. (218) 

Note that these last two relations are possible due to the rotational symmetry of the integrand. Then 
we expand these z-components of momenta upon the hyperspherical harmonic basis using the following 
simple relation (e.g. the p z case): 

Pz = 7—^ — ( ^Y 2 , lfi (n p )) . (219) 



1 + u VV2 

Now we recall the hyperspherical harmonics (see the appendix on hyperspherical harmonics for details) 
that we will be using and their orthogonality and phase relationships: 

y„(fi) = y„,j, m (n) = / n ,,(w)i5, m (M) , (220) 

Y n ,l,m = (-1) +m3 C(,-m ! = ( _ l) m *T*-n» ! /n,J = /nj i (221) 

dfi (4) = = df7 <3) dwsin 2 w , (222) 

J dnY*Y^ = 5^' , y dwsin 2 u>f* J n <j = 8 nn > , j dd ) Y^Yf >m ,i = 5w5 mm > ■ (223) 

After straight-forward application of these relations we obtain: 

,a 4tt r , . 2 (1-cosw) 6tt 2 . , 

• £? = — / dwsin 2 W7; £ = (224 

ei 7 (1 + cosw) ei 

^ 37T 2 ^ 1 f f 71 dus'm 2 uj „ , . „ . .\ 2 ,„„,s 

• ^ = ^E^(/ o (225) 

37r 2 \ ( dwsin 2 w „ , . „ , ,\ 2 ,„„„x 

• ^---E-(y o iT^ /2 ' i(w)/ - i(w) j • (226) 



n=2 

For the Gjj terms, we use the following relations: 

1 1 m 



E-T 11 Bl 2e^ + <^ (22?) 



m m 

^2 



_j_ _ oi*« E ^ (W , M . (228) 
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These give: 



(E-T)G H = ^(i + u ' o ni + u )J2^(n p )Y;(n p/ ) 

1 a 



(229) 



We use the rotational symmetry of the integrand and expand the integrand on the hyperspherical har- 
monic basis as was done for the three Gi terms. Then we have: 



37r 2 ^ 1 / rdu>sin 2 uj 



2ei — ' n 

1 n=2 



3^ 

2ei 



1 / f" dui sin 2 lo 

n 3 



n=2 



ei ^n 2 



1 + COS Ul ' 



du sin 2 w N 2 

-/2,l( W )/n,l( W ) 



1 + COS U! " 



(230) 
(231) 
(232) 



We use the same relations for the Gjjj terms and for the three Gjj terms, and we use the rotational 
symmetry of the integrand to rewrite the appropriate pieces of the integrand in terms of ^2,1,0 as we did 
for the Gi and Gjj terms. However, we need to discuss one additional relation that allows the remaining 
SM 2 (III) calculation to be done analytically. In Schwinger's 1964 paper |Q he gives the following 
formula: 



1 



2tt 2 (1 - pf + p(u - u'f 



OO 1 

-^y nAm (o)y n * i)m (Q'), 



(233) 



n=l 



where u and v! are of unit length and < p < lF^] Inside the brackets in Gm we have: 



1 



1 dp 

7 \Mp(p - P') 2 + C(l - pf C 



1 



1 ff£i 

7 c V(i-p) 2 + p(u-u') 2 



- 1 



(234) 



Recall C = (e 2 + p 2 )(ef + p' ). Also recall that we are using the coordinate change of Eqs. (E03) and 



(203). Eq. (181) with e„ = e\ then applies and was used. In Eq. (234), < p < 1 and u and it' are of 
unit length, thus Schwinger's equation can be used and we have: 

am(l + u ){l + u' a ) 2 i- 1 — « n ~ 2 



(E - T)G IH 



2ej 



dp £ 



P 



YJQ p )Y:(n p ,) . 



(235) 



m#(i,o,o) 



2 This is easily derivable from a more general standard formula that Schwinger gives, 
11 °° P™ 1 1 
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Now, since n > 2 in this sum we can do the integral over p: 



1 n-l 

n-2 _ P 



d PP n 



n- 1 



n - 1 



(236) 



and we obtain: 



(E - T)G in 



I \2 



am(l + u )(l + u' a ) 



•f— [ n(n — 1) p 



^(1.0,0) 



(237) 



For terms in £ which contain Gin, one obtains: 

3^ 2 



_ 3tt 2 ^ i / r 

2ei ^ n(n - 1) V7o 1 + cosu 



dusin 2 w N 2 

/2,l( W )/n,l( W ) 



in 



3^ 2 
2e 



diu sin 2 < 



~~ E "T7~ — TT ( / T "'" — /2,iM/„,i(w) 

i ^ 1 [ doj sin 2 w 
i ^ n 2 (n - 1 



n=2 

3tt 2 

e 



n=2 



1 + COS U) 



/2,l( w )/n,l( W ) 



(238) 

(239) 
(240) 



Now recall \ = fpr£, and also notice that all the summands are the same, thus putting it all together 
we have: 



X 



3 3 
2 + 8 



til 



1 



1 | 1 

T 2 I T 3 qri 3 ( 



1 



— ' In n(n — 1) rt 2 n 3 n 3 (ri— 1) 
n=2 \ v ' v ; 



3 

8 I „ 



2 2 

n n 2 n 2 (n — 1) / \Jq 1 + cosw 



7T T • 2 \ ^ 

sin w 

/2,l( W )/n,lM 



2 8 — \ n 

n— 2 



dio sin 2 ■ 



-/2,i( w )/™a( w ) 



/q 1 + COS U) " 

The remaining sum can be done analytically. To see this, first define two integrals: 



(241) 
(242) 



h 
h 



p p = UttY 

1 + Un 1 + U' 



TT 7 -2 \ ^ 

au> sin oj 



g 1 + cos u> 



16tt 4 



/" dfip dJV (k - k') 2 _ 2 

y i + u 1 + < ( P - p') 2 ~ 3 1 ■ 



(243) 
(244) 



The last equality followed from rotational symmetry of the integrand. We can also calculate I2 a hard 
way which gives :p| 

1 / / aw sin < 



J 2 = 16tt 4 - 4tt 4 V - 

Z — ^ r> 



UJ 

1 + COS Cl> 



/2,lM/™,lM 



32tt 4 



(245) 



We use Eq. (|21S|), = ^ ^ (Y 2 ,i,-i («p) - *2,1,1 («*)) and Py = -j^- f (Y a ,i,-i (Hp) + ^2,1,1 («p)) 
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This last equality followed from Eqs. (|243| ) and fl244|) (the easy way to calculate 12)- Thus, we have 



E 

n=2 



if r* dui sin u> 



, , .. /2,l(w)/n,l( w ) . 

71 \7o 1 + C0S w J o 



Combining this result with Eq. (242) gives 



(246) 



3 3 (A\ 



(247) 



Thus, recalling Eq. (214), we have 



SM 2 



24 



-(351+92) , 



(248) 



where <?i and gi are given in Eqs. (210) and (211) respectively. 



APPENDIX D: M 2 -VS- B„ 



In this appendix we will invert the equation 



Mi 



(2m + B N ) 2 , 



(249) 



and obtain the a-expansion for the binding energy, B N . In this work we set up a procedure to calculate 
M 2 . This gave 

N 



M 2 = M 2 N + b 4 m 2 a 4 + O (a 5 ) 



(250) 



ith 



M 2 N = 4to 2 + 4mS A 



(251) 



For the lowest order spectrum of H R we obtained 



1 ma 2 
4 n 2 



(252) 



Taking a square root of M 2 gives 



iTta 4 / 64 1 



2 32n 4 



(a 5 ) 



(253) 
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Recall that = a 2 Ryd. Now, in this work, Eqs. fll55| )-( |158| ) are the results of our calculation of the 
spin splittings of M 2 in the ground state of positronium. These were derived in the form of Eq. ( p50[ ) 
with result 

7 

b^triplet) — b^singlet) = — . (254) 

3 



Given Eq. fl253| ) this implies 

B triplet -B 3inglet = 7 -a 2 Ryd + 0(a 5 ) . (255) 

This we recognize as the well known result for the positronium system. A final note is that if the physical 
values of the fine structure constant and Rydberg energy ( 13 * Q and 13.60eV respectively) are applied to 
this previous formula, the result agrees with experiment to one-half of a percent ]19[ . 
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Figure captions 

Figure 1: This illustrates the spin and momenta label conventions used in this paper. 

Figure 2: SM^ is the part of the ground state mass spin splittings from Eq. (144). m is the electron 



mass and a is the fine structure constant. The state labels 1, 2, 3 and 4 are explained in Eq. (151). The 
two upper most levels should coincide in a rotationally invariant theory. 

Figure 3: The combined ground state mass spin splitting in positronium to order a 4 is illustrated 
using the same notation as in Figure 2. <5M 2 is given by Eq. (152) and is calculated in Appendix C. The 



final combined result (on the right) corresponds to a rotationally invariant theory. 
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